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ABSTRACT 

We examine continuous measurements of the high-degree acoustic mode frequencies of the Sun covering 
the period from 2001 July to June 2014. These are obtained through the ring-diagram technique 
applied to the full-disk Doppler observations made by the Global Oscillation Network Group (GONG). 
The frequency shifts in the degree range of 180 < £ < 1200 are correlated with different proxies of 
solar activity e.g. 10.7 cm radio flux, the International Sunspot Number and the strength of the local 
magnetic field. In general, a good agreement is found between the shifts and activity indices, and 
the correlation coefficients are found to be comparable with intermediate degree mode frequencies. 
Analyzing the frequency shifts separately for the two cycles, we find that cycle 24 is weaker than cycle 
23. Since the magnetic activity is known to be different in the two hemisphere, for the first time, we 
compute the frequency shifts over the two hemispheres separately and find that the shifts also display 
hemispheric asymmetry; the amplitude of shifts in the northern hemisphere peaked during late 2011, 
more than two years earlier than the south. We further correlate the hemispheric frequency shifts 
with the hemispheric sunspot number and mean magnetic activity index. Since the frequency shifts 
and the hemispheric activity indices are found to be significantly correlated, we suggest that the shifts 
be used as an indicator of hemispheric activity since not many indices are measured over the two 
hemispheres separately. We also investigate the variation at different latitudinal bands and conclude 
that the shifts in active latitudes correlate well with the local magnetic activity index. 
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1. INTRODUCTION 

The eigen-frequencies of solar oscillations which pro¬ 
vide insight into the interior structure and dynamics 
of the Sun are known to vary with the solar activity 


dWoodard fc Noyes|1985 Elsworth et al.||1990 Woodard 

et al. 1991|). With improved observations, the relation¬ 

ship between mode frequencies and activity cycle has 
been extensively studied and it is well established that 
the frequencies exhibit a good correlation with many dif¬ 


ferent solar activity proxies (Chaplin et al. 2007 Jain 


et al. 2009|. Comprehensive analyses using consistent 
and longer data sets both from ground- and space-based 
instruments, however, have shown different and unex¬ 
pected results. For example, a phase-wise analysis of 
the intermediate degree frequencies showed that the de¬ 
gree o f correlation is diff erent from phase to phase of the 
cycle (Jain et al. 2009). The extended minimum phase 
between cycle 23 and cycle 24 also revealed significant 
diffe rences between acousti c frequencies and solar activ- 
itv (iBroomhall et al.||2009[ ISalabert et al.||2009| iTripa- 
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thy et al. 


lations JNe 


2010b r Analyzing Hirmmghami bolar Uscil- 
work data over three solar cycles, |Basu et ah 
(2012) presented evidence that the solar sub-surface layer 


was very different between cycles 22 and 23, with deeper 
changes occurring during cycle 22 compared to cycle 23. 
A recent analysis of low-degree modes from the Global 
Oscillations at Low Frequency instrument also showed 
that cycle 24 is magnetic ally weaker in the u pper sub¬ 
surface layers of the Sun ( Salabert et al.|[2M5 ). 

Since the frequency shifts are predonimantly a func¬ 
tion of the frequency alone, it has been suggested that 
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shallow layers (Elsworth et al. 1994[ 

Howe et al. 

1999 

Chaplin et al. 2001 

Simoniello et al. 

2013b). In addi- 

tion, a small but statistically sigmtican 
with turning points at or below the co 

change in modes 
nvection zone has 

also been reported (Chou & Serebryanskiy 2005 

Bald- 


mams an open question viz. whether the shifts are caused 
by structural changes in the sub-sur face layers, such as 
changes in temperature (|Kuhn|1988|) or by the magnetic 
fields either through the action of the Lorentz force on 
the plasma or indirectly by modifying t he physical prop¬ 
erties e.g. s ize of the acoustic cavity (Dziembowski & 
GoodepOOS ). Anothe r possibility is the inf luence of the 


magnetic atmosphere dJain & Roberts|1996|). The analy- 
f - V,,IIA„l 


sis of frequency shifts by [Dziembowski fe Goodej further 
suggest that it is the indirect effect that dominate, in 
particular changes to the near surface stratification re¬ 
sulting from the suppression of convection by the mag¬ 
netic field. On the other hand, they also reported small 
but significant departures from the simple scaling of the 
lower-frequency p mode frequency shifts with the inverse 
mode inertia implying that there is a contribution from 
deeper layers due to the direct effect of the magnetic 
fields. However, through theoretical modeling, [Foullon 
'& RobertsI (|2005|) concluded that the effect of a buried 
magnetic neld on the frequencies of solar acoustic modes 
is too small to account for t he changes in freque ncy over 


the solar cycle. Moreover, Howe et al. (2002) showed 


that the temporal and latitudinal distribution of the fre¬ 
quency shifts is correlated with the spatial distribution of 
the surface magnetic field implying that the shifts cannot 
be purely explained by the structural changes. Thus, a 


























































































2 


Tripathy, Jain, and Hill 


thorough modeling is required to better understand the 
cause of the frequency shifts. 

In this context, it is desirable to analyze high-degree 
modes whose lower turning points occur very close to 
the solar surface (Figure 0. But determination of high- 
degree mode frequencies is a challenging task since in¬ 
dividual mod es blend into ridges at h igh-degree (see, 
for example, [Korzennik et ah |2013b[ and references 
therein). In spite of this inherent difficulty, there have 
been several attempts to compute and analyze high- 
degree oscillation frequencies. Using time series from 
the Michelson Doppler Imager (MDI) on board the So¬ 


lar an d Heliospheric Observatory {SOHO), Rhodes et al. 
(20021 estimated the high-degree mode frequencies and 


demonstrated the response of the p-modes to short-term 


changes in solar activity. Rose et al. (20031 extended 


this analysis to include data from Global Uscillation Net¬ 
work Group (GONG). Considering global modes from 


MDI, Rabello-Soares et al. (2008b) concluded that the 
changes in the high-degree mode frequencies are in good 
agreement with the intermediate-degree results, except 
for years when the MDI was highly defocussed. A sim¬ 
ilar study, but including more data, demonstrated a 


tivitv indices (iRabello-Soares 

2011 

, which was not con- 

firmed by 

I'ripathy et al. ( 

20 

[3b). 

Rhodes et al. 

(2011) 

also investigated the tempora 

changes in high-degree p- 


mode frequencies and pointed out that short-duration, 
high-frequency shifts are more sensitive to changes in so¬ 
lar activity than the shifts measured by long-duration 
studies consisting only the lower frequencies. Recently, 


Korzennik et al. (2013a) computed high-degree modes for 
three ditterent epochs from three different instruments 
and compared the mode parameter s. 


Ring-diagram analysis (Hill|1988), in addition to global 


helioseismic techniques, measures the characteristics of 
high-degree modes (180 < £ < 1200) in localized regions 


around active regions (IRaiaguru et al. 12001 iBogart et al. 

2002 Rabello-boares et al.|20U8a). 

Hindman et al. 
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netic activity which was lat er shown to be less reliable 
during low activity periods (Tripathy et al.|[2MT). How¬ 


ever, the response of the frequencies ot the high-degree 
modes to solar activity using a continuous and consistent 
data set has not been undertaken previously. The results 
reported earlier mostly use discrete data sets obtained by 
the MDI instrument when it was operating in its high- 
resolution mode three months every year. Our goal, in 
this paper, is to investigate the behavior of high degree 
mode frequencies over a period of about 13 years con¬ 
sisting of both the descending phase of cycle 23 and the 
ascending phase of cycle 24. This analysis supplements 


our previous study (Tripathy et al. 2013a hereafter pa¬ 
per I) where we had primarily focused on the minimum 
phase of solar cycle 23. Since magnetic activity (e.g. 
the sunspot number and polar flux) exhibit north-south 
hemispheric asymmetry, for the first time, we examine 
the variations in oscillation frequencies over these two 
hemispheres separately. We further explore the relation¬ 
ship between the frequency shifts and magnetic activity 
at different latitudes in both hemispheres. 

The rest of the paper is organized as follows. The data 
and the analysis procedure is described in § 2, while the 


results are described in § 3. Finally, we summarize the 
main conclusions in § 4. 

2. DATA ANALYSIS 

The Global Oscillation Network Group has been op¬ 
erating in high-resolution mode since July 2001 and 
produces full disk Doppler images of the Sun every 
minute. These images are then processed through the 
ring-diagram pipeline (Gorbard et al. 2003) to produce 
the mode frequency tables. In short, the ring-diagram 
technique ( Hilljl988 ) produces a three dimensional spec¬ 
tra (two in space, one in time) of a small region on the 
surface of the Sun tracked and remapped over a given 
time period. In the GONG pipeline, each region (here¬ 
after tile) covers an area of 15° x 15° in heliographic lon¬ 
gitude and latitude and is tracked for a period of 1664 
minutes (hereafter ring-day). The centers of the tiles 
are separated by 7.5° and a set of 189 overlapping tiles, 
referred to as dense-pack, covers the region of ±52.5° 
in latitude and central meridian distance from the disk 
center. The three di mensional spectra i s fitted with a 
Lorentz profile model (Haber et al.||2000| 


P{kx,ky,v) = 


A 


■ + TA> (1) 


(z.-i/o±A:.C/,±fcyl/j,)2+F2 fc3’ 

where P is the oscillation power for a wave with a tempo¬ 
ral frequency ly and the total wave number k'^ = k^ + ky. 
There are six parameters to be fitted: the zonal and 
meridional velocities, 14 and Uy respectively, the cen¬ 
tral frequency, ly, the amplitude. A, the mode width, F 
and the background power, b. For each tile, the peak fit¬ 
ting procedure returns about 200 modes with degree, £, 
between 180 and 1200 and with frequency between 1700 
and 5600 /rHz. It is, however, to be noted that these are 
not discrete modes in the same sense as those observed 
at lower degrees. The ring frequencies are taken from the 
GONG websit^and cover the period from 2001 July 26 
to 2014 June 19. 

Since GONG is a ground-based network of six stations, 
poor observing conditions and instrument failures give 
rise to gaps in the data. These gaps are measured in 
terms of a fill factor, where a fill of unity means no data 
loss (100% coverage). The distribution of the fill factor 
of the ring data are shown in Figure Since gaps in 
the data causes biases and data of lower fill may not be 
fitted reliably, we exclude data where the fill factor is 
lower than 70%. This resulted in 3761 out of a total 
of 4146 ring-days and consists of approximately 710,800 
tiles . 

The mode frequencies (and other mode parameters) 
measured by the ring-diagram technique are affected by 
foreshortening as well as gaps in obse r vation . Following 
the method outlined by Howe et al. (2004), we model 
the effects of position of the disk as a two-dimensional 
function of the distance from disk center p without the 
cross-terms, combined with a linear dependence on the 
fill factor, f{t), 


v { Px : Py : t ) = Oq ±01/9^; ±02/3^ TasPy ±04^^ ±05/(1), (2) 

^ http://nisp.nso.edu 







































































Acoustic Frequencies of the Sun during Solar Cycle 23 and 24 


3 



0 200 400 600 800 1000 1200 

Spherical Harmonic Degree (2) 



Fill {%) 


Figure 1. The lower turning point of oscillation modes for inter¬ 
mediate (blue circles) and high (red circles) as a function of the 
spherical harmonic degree, £. The spread is due to different radial 
orders. 


Figure 2. Fill factor of GONG ring-data between 2001 July 26 
and 2014 June 19. 


where px and py are the longitudinal and latitudinal com¬ 
ponent of p and ai are coefficients determined by fitting. 
In order to minimize the effect of the magnetic field on 
the fitting, we use the frequencies over a two year period 
when solar activity is low (between 1 September 2007 
and 31 August 2009) and evaluated the a^s for each (n, 
i) mode. The resultant coefficients and the resultant dif¬ 
ferences between the corrected and original frequencies 
are shown in Figures 3 and 4 of paper I. 

Figure compares the fitted modes obtained from 
global and ring-diagram analysis techniques in a char¬ 
acteristic i — v plane. The black circles display the inter¬ 
mediate degree modes fitted using MDI data corre spond¬ 


ing to the epoch 2002 March 31 to 2002 June 10 (Schou 
1999|). The red circles represent the global high-degree 


modes from the GONG data obtainecQfrom a time series 


averaged over the same 98 days; the lowest ridge corre¬ 
sponds to the / mode while others are the p-mode ridges 
for radial order n from 1 to 6. It is evident that the fre¬ 
quencies computed from different analysis mostly agree 
with each other, as the symbols appear to fall along the 
same ridges. This also suggests that the systematic er¬ 
rors in the frequencies arising from analysis choices are 
independent of time. 

We further compare the individual frequencies ob¬ 
tained from GONG data between high- (2002) and low- 
(2010) activity periods where the frequencies correspond¬ 
ing to low activity period are computed from a 67 day 

^ https: //www. cf a. harvard. edu/$\backsiin$sylvain/ 

research/tables/HiL 



14 ( 

Korzennik et al. 

2013a 

). The blue circles display the 
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high-degree modes o 

stained by the ring-diagram analysis 

Spherical Harmonic Degree 


Figure 3. The l—u diagram of the modes obtained from the use 
of different techniques and data source. The black circles display 
the fitted intermediate degree modes using MDI data, the red dots 
represent the global high-degree modes from GONG, and the blue 
circles display the high-degree modes obtained through the ring- 
diagram analysis of the GONG data. 


time s eries star^g from 2010 May 7 (Korzennik et al. 
2013a I. Figure |4| shows the frequency differences. The 
pluses denote the differences between the global high- 
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degree modes while the filled circles shows the differences 
between the modes obtained from the ring-diagram anal¬ 
ysis. The mode frequencies obtained from two different 
methods agree very well except for the scatter seen in 
the global modes at high frequencies. 



1000 2000 3000 4000 5000 6000 

Frequency (/uHz) 


Figure 4. Frequency differences between 2002 and 2010. The plus 
symbols denote the differences computed from high-degree global 
modes obtained from Korzennik’s fitting in the degree range of 100 
and 1000. The filled circles denote the differences between high 
degree modes computed using the ring-diagram technique. The 
2002 frequencies are calculated from a time series of 98 days with 
a start day of 2002 February 14, while those in 2010 are computed 
from a time series of 67 days starting from 2010 May 7. 


For the purpose of characterizing activity associated 
with each tile, we compute a magnetic activity index 
(MAI) which indicates the local magnetic field strength 
associated with each tile. The index is calculated by 
tracking and remapping the line-of-sight magnetograms 
during the same interval as the Dopplergrams and av¬ 
eraging the absolute values of all the pixels with a field 
strength higher than a threshold value. Since the entire 
period under study is not covered by the MDI magne¬ 
tograms alone, we use the 96 minute MDI magnetograms 
up to 2007 April 9 and GONG magnetograms sampled 
every 32 minutes thereafter. A lo wer threshold val ue of 
50 G and 8.8 G are used for MDI (Basu et al.|2004 1 and 
GONG respectively, which are approximately 2.5 times 
the estimated noise in individual magnetograms. For 
generating a uniform magnetic activity index, the MDI 
MAI values are scaled to GONG values using the con¬ 
version factor of 0.31 that was obtained by relating the 
line-of-sight magnetic field measurements between these 
two instruments u sing a histogram-equalizing technique 
( Rfiey et al. [2014 1. 

A well known systematic effect in solar physics and 
specially with ground-based resolved data is the annual 
variation with the solar inclination toward Earth which 



Year 


Figure 5. Temporal evolution of the mean magnetic activity in¬ 
dex, MMAI. The solid black line shows the running mean over 181 
points. The dashed red line represents the 181 day smoothed 10.7 
cm radio flux measurements scaled linearly with MMAI. 


can lead to a loss in spatial resolution and hence change 
the geometric foreshortening. This effect is clearly seen 
in the MAI, the frequency shifts (calculated later) as well 
as in solar activity indices. We have applied a technique 
based on Fourier Transform to remove the 1 year peri¬ 
odicity from all the date sets used in the study of tem¬ 
poral variations. However, for a better demonstration, 
uncorrected data has been included in some figures by 
means of symbols. Also note that the frequency shifts 
and MAIs have been interpolated to 4146 ring-days be¬ 
fore the Fourier Transform was applied. 

Figure shows the temporal variation of the mean 
MAIs (MMAI) where the mean is taken over 189 dense- 
pack tiles of each ring-day. We also calculate a 181 point 
running mean, which corresponds approximately to seven 
Carrington rotations in terms of ring-days, to smooth 
out the short-term variations. This is shown by solid (or 
dashed) lines in most of the figures. In the same plot, 
we also show the variation of the 10.7 cm radio flux, 
Eioy as a measure of the solar activity (dashed red line). 
Significant agreement is seen between the shape of the 
curves indicating a high correlation which is confirmed 
with a rank correlation coefficient value of 0.97. Thus 
the MMAIs will be considered as an activity proxy in 
rest of the analysis. 

In order to examine the sensitivity of the frequency 
shifts to solar activity, we further use two commonly 
used activity proxies; the International Sunspot Num¬ 
ber, ISN, and Fio.y both obtained from National Geo¬ 
physical Data Center (NGDC)R The sunspot number is a 
weighted estimate of the individual sunspots and sunspot 

® http://www.ngdc.noaa.gov/stp/solar/solardataservices. 
html 
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groups derived from inspection of the solar photosphere 
in white light while the radio flux index is a measurement 
of the integrated emission at 10.7 cm wavelength from all 
sources present on the solar disk. 

3. RESULTS AND DISCUSSIONS 
3.1. Spatial Variation 


of 


Hindman et al. (2000) studied the spatial variation 
the frequencies obtained from ring-diagram analysis 
and showed that the active regions appear as locations 
of large frequency shifts. This gave rise to the conclu¬ 
sion that the frequency shifts can be used as a tracer of 
surface magnetic a ctivity. Subsequent analysis (Tripathy 
et al. 2010a 2011), however, indicated that the conclu- 
sion IS valid only during the period of high activity and 
breaks down during the activity minimum period. 

We estimate the spatial frequency shift, 5vs, weighted 
by the fitting error, for each tile according to 


Siy,{t) = 


(3) 


where the frequency difference of each mode, is 

calculated with respect to the error weighted average of 
that particular mode present in the dense-pack of each 
ring-day. A sequence of images showing 5vs and MAI for 
four days during the high- and l ow-activity periods are 
presented in Figures 1 and 2 of Tripathy et al. (2011) 
demonstrating the relationship between them. 

For a quantitative analysis, we calculate the Pearson’s 
linear correlation coefficient, rp, between Si's and MAI 
of each tile in the dense-pack. The resultant coefficients 
(color coded for different years) as a function of the 
MMAIs are shown in Figure |^. We observe that for 
high values of MMAI the correlation is stronger while 
for MMAI < 1 G, the correlation coefficients are small 
with a median value of 0.41. On closer inspection, we 
also find higher correlation coefficients for small MMAI’s 
specially during the minimum phase. For a better visu¬ 
alization, we plot the correlation coefficients for different 
phases separately. Figure]^ shows the period around the 
maximum phase of cycles 23 (2001-2003) and 24 (2012- 
2014) and illustrates stronger correlation between the fre¬ 
quency shifts and mean magnetic activity index. The 
coefficients for the descending phase of cycle 23 (2004- 
2006), the extended minimum phase (2007-2009), and 
the ascending phase of cycle 24 ^010-2011) are shown 
in panels c, d and e of Figure |m respectively. In all 
cases, we note both high and low correlation coefficients 
for MMAI values of « 1 G and a large spread in the val¬ 
ues of the coefficients (largest for the extended minimum 
phase). We further observe negative coefficients in pan¬ 
els c and d (fewer in panel d) implying anti-correlation 
but the values are not significant. Analysis of global 
mode frequencies with solar activity proxies during the 
extended minimum phase between cycles 23 a nd 24 has 
also indicated anti-correlation bet ween them (Salabert 
et al.||2009| [Tripathy et al.||2010bj). Thus our result is 
consistent with the global low- and intermediate- degree 
modes. 

The distribution of the linear correlation coefficient in 
the form of a histogram is shown in Figure Only a 
small fraction of 63 ring-days (1.7%) are anti-correlated 


with MMAI while 73% and 65% ring-days have positive 
correlation coefficients higher than 70% and 80%, respec¬ 
tively. 

3.2. Variation with Solar Cycle 

The frequency variation with the solar cycle is investi¬ 
gated by calculating a mean frequency shift, i5u, using a 
formula which is generally applied to compute frequency 
shifts of global modes, 




El=l SVnAt) 


1 




/E;A. 


(4) 


C,e 


where 6iyn,e is the frequency difference of each (n, £) 
mode with respect to the average frequency of the same 
mode over the dense-pack tiles corresponding to a quiet 
ring-day 2011 May 2008. Table H shows the values at¬ 
tained by the mean frequency slmts and three activity 
proxies during cycles 23 (up to 2008 December 31) and 
24. However, the values presented here may not be the 
true minimum and maximum values of the solar cycles 
(and hence the ratio) due to the period considered in 
this study. The final column in the table shows the ra¬ 
tio of amplitudes of cycle 23 to cycle 24, in the sense 
(max-min) 23 /(max-min) 24 . While the mean shifts and 
activity indices show a higher amplitude in cycle 23, the 
MMAI, which is the strong field component of the mag¬ 
netic flux, shows a higher amplitude in cycle 24. This 
is also evident in Figure where Sn and activity prox¬ 
ies are plotted as a function of time. The solid line in 
each panel represents the running mean over 181 points 
and guides the eye to visualize the descending phase of 
cycle 23, ascending phase of cycle 24 and the extended 
minimum period between the two cycles. While the ex¬ 
tended minimum phase is unambiguously detectable in 
the activity indices, it is not noticeable in the frequency 
shifts illustrating that (i) the change in frequency with 
activity is no longer linear during the solar minimum 
and (ii) the minimum phase in frequency lasted over a 
shorter period compared to the solar activity. This is 
consistent with our earlier conjecture that the changes in 
high-degree frequencies may be associated with two com¬ 
ponents: strongly localized acti ve regions and the smal l 
scale horizontal magnetic field (Tripathy et al. 2010a). 
The former explains the larger shifts, while the later is 
responsible for the smaller shifts in all phases of the solar 
activity includ ing t he minimum phase. We return to this 
discussion in § |374| 

We further examine the sensitivity of the frequency 
shifts to solar activity by performing a linear fit be¬ 
tween them using data over the entire 13 year period 
as well over different phases separately (Figure]^. The 
straight lines are the result of linear regressions for cy¬ 
cle 23 (blue), cycle 24 (green) and all data (red). The 
black solid line in each of the panels is the result of a 
quadratic fit using all data and differs from the linear 
fit (red line) when activity is high. Thus, we conclude 
that the activity proxies display simple linear relations 
with the frequency s hifts instead of a qua dratic relation 
as reported earlier Rabello-Soares 11 2011 1) . We conjec¬ 


ture that the quadratic relation is an outcome of using 
fewer data points as we observe at high-activity period. 
Table shows the frequency shift per unit activity gradi¬ 
ents and the corresponding linear correlation coefficients 
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Figure 6. (a) Pearson’s linear correlation coefficient, Vp, between 5us and MAI over 189 dense-pack tiles for each ring-day as a function 
of the mean MAI. The symbols and colors represent the year-wise variation and are indicated in the top panel. Panels (b-e) shows 
the coefficients around the maximum (2001-2003; 2012-2014), descending (2004-2006), minimum (2007-2009) and ascending (2010-2011) 
phases of the solar activity in cycles 23 and 24, respectively. Note that the range of Mean MAI in each panel is different. 
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Figure 7. Histogram showing the distribution of the correlation 
coefficients as a function of the number of days. 


for each of the proxies. For MMAI, the gradient is higher 
for cycle 23 compared to cycle 24, while for other proxies 
the gradients do not vary significantly between different 
cycles or over the entire time span. The level of cor¬ 
relation with different activity proxies are higher than 
80% largely in good agreement with the analysis of in¬ 
termediate degree global modes ( Jain et ^|2011 1; there 


are, however, marginal changes between cycle 23 and 24 
and entire period. Also, as anticipated, the correlation 
is higher between 5v and magnetic activity measured lo¬ 
cally. 

Before the data were corrected for the annual varia¬ 
tion, a secondary periodicity except during the minimum 
phase was observed in the frequency shifts as well as in 


biennial signal 1 

Broomhall et al. 

2012 

Simoniello et al. 

2013b 

) that has been recently receivec 

a great deal of 


3.3. Hemispheric and Latitudinal Variation 

In addition to the temporal changes, magnetic activity 
is generally north-south asymmetric and varies with lat¬ 
itudes. The hemispheric asymmetry in sunspot number 
and polar flux is well established and provides important 
clues for und erstanding the mech anism of solar cycle (for 
a review see, INorton et al. 2014). The variation can be 


easily seen in magnetic butterfly diagrams where the ac¬ 
tivity belts move towards the geometric equator with the 
progression of the solar cycle. In addition, the wings of 
the butterfly diagram show a high degree of mirroring 
across the equator indicating coupling between the two 
hemispheres. In this context, the frequencies obtained 
from ring-diagram technique provides a natural way to 
investigate variations between different hemispheres as 
well as latitudes since the m-averaged global frequencies 


of oscillation modes, commonly used for solar cycle stud¬ 
ies, are latitudinally averaged and can not distinguish be¬ 
tween the northern and southern hemispheres. Nonethe¬ 
less, the frequency of individual multiplets (which are 
function of n, £, m) can be be used t o study latitudinal 


variations based on m/£ values (e.g. Howe et al. 2002 


Jain et al. 2011). The sectoral and near-sectoral modes 


are sensitive to the region near the equator while the 
zonal and tesseral modes are sensitive to the wider range 
of latitudes. 

Figuref^shows the mean frequency shifts for northern 
(blue circles) and southern (orange circles) hemispheres 
and the solid line represents the running mean over 181 
points. The shifts are calculated using Equation (3) 
where di'n,i corresponds to the frequency difference of 
each mode with respect to an average frequency over the 
same hemisphere. Both the northern and southern fre¬ 
quency shifts decreased starting from mid-2001 and their 
variations are similar until the end of 2007. Around the 
beginning of the cycle 24, the two shifts diverged and 
reached their maximum amplitude at two different peri¬ 
ods. In cycle 24, the northern shifts peaked in late 2011, 
more than two years earlier than the south, as indicated 
by the dashed and dot-dashed lines, respectively. This 
closely agrees with the hemispheric asymmetry seen in 
other activity indices e.g. the MMAI (panel b) and the 
hemispheric sunspot number (HSN, panel c), where the 
latter is obtained from the NGDC. Conspicuously, both 
the dashed and dot-dashed lines are aligned between dif¬ 
ferent panels, illustrating that the frequency shifts very 
closely follow the changes in activity indices over both 
hemispheres. To quantify the relationship, we compute 
the correlation coefficient between them and the shift per 
unit activity gradients. These are tabulated in Table ^ 
The relative change in the gradients points out that the 
shifts remained approximately constant between the two 
hemispheres. Moreover, the Pearson’s linear and Spear¬ 
man’s rank correlation coefficients further confirms the 
linear relation between the shifts and activity indices. 
The values mentioned in the brackets are the coefficients 
corresponding to the smoothed values and, as expected, 
signify a closer interaction. In the absence of many activ¬ 
ity proxies for different hemispheres, we suggest that the 
frequency shifts may be used as an hemispheric activity 
proxy. _ 

Figure shows the mean frequency shifts for a few se¬ 
lected latitudes where the shifts are obtained by subtract¬ 
ing the average frequency computed at the same latitude 
covering all data sets but only over the central meridian. 
For clarity, we only show the mean shifts resulting from 
a running mean over 181 points. The shifts show sig¬ 
nificant variations between different latitudes as well as 
between the two cycles. Moreover, a clear extended min¬ 
imum phase is not seen at any of the latitudes. In cycle 
23, peak amplitudes are observed around the beginning 
of 2002 for latitudes 7.5° N and 15° N, while for south¬ 
ern hemisphere the peak amplitude can be seen for all the 
latitudes around the same time. On the contrary, we no¬ 
tice opposite behavior in cycle 24; the peak amplitude is 
seen around the end of 2011 for all latitudes while in the 
southern hemisphere, the shifts peak around the begin¬ 
ning of 2014 except for 30° latitude. At this latitude, a 
peak amplitude is observed around 2012 November and 
a decreasing trend afterwords. Interestingly, the shifts 
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Year 


Figure 8. Temporal evolution of (a) 5u, (b) MMAI, (c) F 10.7 and (d) ISN. The solid black line in each panel represents the running mean 
over 181 points. The error bar on the top left side of panel (a) represents 20 times the mean error. 
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Figure 9. Correlations between the frequency shifts and different activity indicators for cycle 23 descending (filled circles) and cycle 24 
ascending (open circles) phases. The straight lines are the linear regressions for cycle 23 (blue), cycle 24 (green) and all data (red). The 
activity proxies are given on the abscissae. The black line is the result of a quadratic fit using all data. 
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Figure 10. Temporal variation of (a) Su, (b) magnetic activity and (c) hemispheric sunspot number (HSN) for northern hemisphere (blue 
circles) and southern hemisphere (orange circles). The solid line in each panel represents the running mean over 181 points, the blue and 
red lines denoting northern and southern hemispheres, respectively. The dashed and dot-dashed lines indicated the peak amplitudes in 
northern and southern hemispheres, respectively. 
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Table 1 

Values of Frequency Shifts and Activity Proxies for Different Solar Cycles 


Activity Index 

Cycle 23 

Cycle 24 

(Max — Min) 23/(Max — Min) 24 

Minimum 

Maximum 

Minimum 

Maximum 

Sv 

- 0.52 

2.24 

- 0.42 

1.99 

1.14 

MMAI 

0.23 

5.80 

0.28 

7.01 

0.83 

Fi0.7 

76.2 

281 

80.4 

253 

1.19 

ISN 

0 

200 

0 

178 

1.12 


Note. — The MMAI has units of Gauss; the 10.7 cm radio flux is expressed in sfu (1 sfu =10 J s ^ 
Hz~^); ISN is dimensionless and Siy' is in /iHz. 


Table 2 

Gradients and Linear Correlation Coefficients of Frequency Shifts as a Function of Activity Proxies 


Activity Index 

Cycle 23 


Cycle 24 


All Data 


Gradient 

rp 

Gradient 

rp 

Gradient 

rp 

MMAI 

0.355 ± 1.58x10-* 

0.91 

0.242 ± 1.60x10“* 

0.90 

0.289 ± 1.09x10“! 

0.88 

FiO.7 

0.009 ± 3.99x10“® 

0.87 

0.005 ± 3.14x10“® 

0.79 

0.009 ± 3.45x10“® 

0.87 

ISN 

0.009 ± 4.22x10“® 

0.87 

0.009 ± 4.89x10“® 

0.79 

0.009 ± 3.45x10“® 

0.87 


Note. — Units of the gradients are /iHz per unit activity, where the activity units are defined in Table 


Table 3 

Values of Activity Proxies, Frequency Shifts and Correlation Statistics between Si^ 
and Activity Indices (Pearson Linear Correlation, Vp, and Spearman rank 
Correlation, r^). 


Activity Index 

Min 

Max 

Gradient 

rp 

rs’^ 

MMAI (G) 

0.04 

8.09 

Northern Hemisphere 

0.285 ± 1.61x10“* 

0.90 (0.97)® 

0.90 (0.97) 

HSN 

0 

110 

0.022 ± 1.32x10“® 

0.84 (0.98) 

0.84 (0.97) 

MMAI (G) 

0.00 

0.00 

Southern Hemisphere 

0.284 ± 1.34x10“* 

0.86 (0.94) 

0.87 (0.95) 

HSN 

0.0 

126 

0.022 ± 1.09x10“® 

0.86 (0.98) 

0.85 (0.98) 


^ The calculated probabilities of having null correlations are 0. 

^ The values in the brackets are the correlation coefficients corresponding to the 181 point 
running mean. 
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Figure 11. Temporal variation of frequency shifts as a running mean over 181 points at selected latitudes on northern (panel a) and 
southern (panel c) hemispheres. Line styles have the following meaning: red (dash-dot-dot-dot): 7.5°, green (dotted): 15°, and blue 
(dashed): 30°. Panel (b) shows the shifts corresponding to the equatorial tiles. 
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Figure 12. Variation of linear regression slopes (circles) and Pear¬ 
son’s correlation coefficients (squares) for different latitudes. 


corresponding to the equatorial belt has only one strong 
peak amplitude a 2002 and a highly reduced amplitude 
(about half of cycle 23) for cycle 24. Thus, if we consider 
frequency shifts that correspond to equator we will con¬ 
clude that cycle 24 is weaker compared to cycle 23 while 
at other latitudes, the shifts are nearly similar (slightly 
higher for southern hemisphere) between cycles 23 and 
24. 

The frequency shift per unit magnetic activity (filled 
circles) and the Pearson’s linear correlation coefficients 
(filled squares) between MAI and Siy for different lati¬ 
tudes are shown in Figure Higher values of the gra¬ 
dients are found corresponding to the solar active lati¬ 
tude belt i.e. about ± 30°. In the same latitudes, mag¬ 
netic activity as measured by MAI are also found to be 
larger and hence a strong correlation is observed between 
the MAI and Sty. The correlation coefficient decreases to 
lower than 0.8 at latitudes higher than ± 37.5° where the 
presence of the active region is sparse. From the nature 
of the two curves, it appears that the linear regression 
slopes and the magnitude of the correlation are well re¬ 
lated. A similar result has been reported fo r frequency 
shifts in global modes (Tripathy et al. 2007). For a ro¬ 
bust comparison, it will be worth analyzing the entire 
solar cycle 24 in terms of the ascending and descending 
phases when the data is available. 


3.4. Minimum Phase 

In order to better characterize and quantify the differ¬ 
ences between the frequency shifts and activity proxies 
during the minimum phase, we calculate the length of 
the minimum period by defining a quiet phase as: 

{maximum(x) — minimum{x)) * 0.07 + minimum{x), 


where x represents either frequency shifts or activity 
proxies. The factor of 0.07 is chosen based on the ob¬ 
served quiet phase in ISN and F 10 . 7 . The temporal evo¬ 
lution of <5u, MMAI, Fio ,7 and ISN showing the extent 
of the minimum phase by filled green color is illustrated 
in Figure |13| It is evident that the minimum phase in 
Sv lasted for a brief period compared to the activity. 
The minimum and maximum values of these parameters 
along with the start time, end time and the length of the 
quiet phase are listed in Table The minimum phase in 
frequency shifts lasted for a period of about 16 months 
while all the activity proxies lasted for a period of about 
30 months or longer. We have also computed these pa¬ 
rameters separately for the northern and southern hemi¬ 
spheres. The values are listed in Table Conhrming 
to known results, the activity (and hence the shifts) are 
found to be stronger in the southern hemisphere as com¬ 
pared to the northern hemisphere. We also note that 
the minimum phase lasted longer in northern hemisphere 
compared to the southern hemisphere, probably, an indi¬ 
rect manifestation of the stronger activity in the southern 
hemisphere. We also find an exception in the northern 
hemisphere where the length of the quiet phase seen in 
5v lasted slightly longer as compared to the magnetic 
activity. 

We further computed the length of the quiet phase for 
each latitude. The starting and ending epoch and the 
length of the minimum phase for 5v and MAI are listed 
in Table values above ± 37.5° latitudes are omitted 
since reli^le measurements could not be obtained prob¬ 
ably due to small variations between the minimum and 
maximum phase. Confirming to the earlier findings, the 
span of the quiet phase seen in frequency shifts is found 
to be reduced as compared to the magnetic activity. 

4. CONCLUSIONS 

In this work we have analyzed high-degree mode fre¬ 
quencies of the Sun over a period spanning about 13 years 
to understand their behavior with the solar cycle. These 
frequencies are obtained from the application of ring- 
diagram technique to the high-resolution Doppler obser¬ 
vations made by the Global Oscillation Network Group 
and covers the period of 2001 July 26 to 2014 June 19. 

We compared the mean frequency shifts with three ac¬ 
tivity proxies- the 10.7 cm radio flux, the International 
Sunspot Number and the strength of the local magnetic 
field measured by the magnetic activity index. In gen¬ 
eral, a good agreement is found between the shifts and ac¬ 
tivity indices and the correlation coefficients are found to 
be comparable with intermediate degree mode frequen¬ 
cies. Although we do not have data covering all of cycle 
23 or cycle 24, the phase-wise variations of the shifts in¬ 
dicated that the frequencies responded to solar activity 
equally well in both the cycles although the shifts per 
unit activity hinted that the cycle 24 is weaker than cy¬ 
cle 23. We also noticed that the frequency shifts do not 
show a prolonged minimum phase as is seen in many 
activity indicators suggesting weak correlations between 
the shifts and activity indices during the activity mini¬ 
mum period. 

For the first time, we also computed frequency shifts 
over the north and south hemispheres separately and an¬ 
alyzed them as a function of hemispheric sunspot number 
and magnetic activity index. The asymmetric nature of 
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Figure 13. Temporal evolution of (a) Su, (b) MMAI, (c) Fio.7 and (d) ISN showing the extend of the solar minimum phase (filled area) 
which are overlaid on the running mean over 181 points. 
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Table 4 

Minimum and Maximum Values of Different Parameters, Start Time, 
End Time and Length of the Quiet Phase. 


Index 

Min 

Max 

Start 

End 

Length 




All Data 



Su 

100.2 

1455.4 

2007 Jun 

2008 Oct 

17 

MAI 

0.8 

5.0 

2007 May 

2009 Oct 

30 

FiO.7 

82.8 

241.7 

2007 Feb 

2010 Mar 

38 

ISN 

0.0 

159.5 

2007 Apr 

2010 Jan 

34 



Northern Hemisphere 


Su 

180.8 

1324.4 

2007 Jul 

2010 Sep 

39 

MAI 

0.9 

4.8 

2006 Aug 

2009 Jun 

35 

HSN 

0.0 

80.5 

2006 Mar 

2009 Sep 

43 



Southern Hemisphere 


Su 

28.7 

1700.0 

2007 Jun 

2008 Oct 

17 

MAI 

0.7 

5.6 

2008 Mar 

2010 Feb 

24 

HSN 

0.0 

92.9 

2007 Sep 

2010 Jun 

34 


Note. — Units of Su is nHz, MAI is Gauss and E10.7 is sfu. The length 
is measured in month. 


Table 5 

Start Time, End Time and Length of the Quiet Phase at Various Latitudes. 


Latitude 


MAI 



5y 


Start 

End 

Length 

Start 

End 

Length^ 

0.0° 

2008 Jun 

2011 Apr 

35 

2009 Mar 

2011 Dec 

34 

7.5° 

2006 Apr 

2010 Jun 

51 

2006 Dec 

2010 Sep 

48 

15.0° 

2006 Feb 

2009 Sep 

44 

2006 Apr 

2009 Sep 

42 

22.5° 

2006 Jan 

2009 Jan 

37 

2006 Feb 

2008 Jul 

30 

30.0° 

2006 Mar 

2008 Sep 

31 

2005 Oct 

2008 Mar 

30 

37.5° 

2007 Jan 

2009 Sep 

33 

2005 Aug 

2008 Jan 

30 

-7.5° 

2008 Aug 

2011 Jun 

35 

2009 Jan 

2011 Sep 

33 

-15.0° 

2008 Mar 

2010 Jul 

29 

2008 Mar 

2010 Jul 

29 

-22.5° 

2006 Apr 

2009 Oct 

41 

2006 Oct 

2010 Jan 

40 

-30.0° 

2006 Apr 

2009 Apr 

37 

2006 Jan 

2008 Dec 

36 

-37.5° 

2007 Mar 

2009 May 

27 

2007 Jan 

2008 Dec 

24 


^ Length is measured in month. 
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the magnetic activity is also reflected in the frequency 
shifts. We note that the shifts behaved similarly in the 
northern and southern hemisphere until the end of 2007 
and diverged thereafter. This indicates that the Sun had 
started showing deviations even before the onset of the 
minimum phase and that the mode frequencies sensed 
these changes before it was manifested otherwise at the 
visible surface. Further, we find that the shifts reached 
their maximum amplitude at two different periods-the 
amplitude of shifts in northern hemisphere peaked dur¬ 
ing late 2011 more than two years earlier than the south. 
Such asymmetric behavior between the ascending and 
descending phase of sola r cycle 23 has been studied by 


Simoniello et al. (2013a) and it is argued that an in- 
crease of asymmetry might be used as a precursor of 
a grand minima. Theoretical calculations also indicate 
that strong asymmetries are induced by the quadrupolar 
component and characterize per iods before entering or 
leaving the grand minima state (Tobias|[l996 Knobloch 


et al.|19^ ). Thus in our view the asymmetry associated 


with the frequency shifts may have certain implications 
for the dynamo action taking place either in the sub¬ 
surface shear layer or at the base of the convection zone. 

Since the shifts and the hemispheric activity indices 
are found to be significantly correlated, we suggest that 
the frequency shifts may be used as an indicator of hemi¬ 
spheric activity since not many indices are measured over 
the two hemispheres separately. We also investigate the 
variations at different latitudes and conclude that the 
shifts in active latitudes correlate well with the local 
magnetic activity index. We believe that this analysis 
provides impetus to analyze solar cycle 24 in terms of 
hemispheric asymmetry and different phases to compre¬ 
hend the magnetic cycle of the Sun. 
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